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ABSTRACT 



The determination of the effective mass and the dairying of 
marine propellers is essential to the analysis of the longitudinal 
vibrational characteristics of a ship’s shafting system. Effective 
mass is the actual mass of the propeller plus a virtual mass due to 
entrained water vibrating with the propeller. Our object was to con- 
tribute to the solution of the general problem by making an experiment- 
al determination of the virtual mass under operating conditions. 

Damping was to be measured coincidently. 

Propellers selected for testing included a set of merchant type 
propellers with P/D values ranging from 0.8 to 1,2 (U bladed), a set 
of DD type propellers with a variation of P/D (1* bladed), and a set of 
DD type propellers with a variation of P/D (3 bladed). Operating condi- 
tions (in the MIT Propeller Tunnel) were attained by driving the propel- 
lers through a torque coupling having small longitudinal stiffness and 
maintaining appropriate water velocities. The propellers were supported 
by a shaft and housing which were principally supported by two sets of 
pairs of converging rods, so arranged as to allow longitudinal motion 
of the housing but constraint against any other movement. A step longi- 
tudinal displacement of the propeller and housing caused a transient 
vibration which was recorded using a strain-gage transducer and recorder. 
The primary longitudinal stiffness of the system was related to the 
stiffness of the rods in bending. A longitudinal stiffness constant was 
determined which related the mass of the whole system to the square of 
the natural period. Using this relationship the mass could be deter- 
mined by observing the period of vibration. Th« original decay was 
measured from the tape record of the transient vibrations. 

The virtual mass was determined for each propeller along with 
the damping. The data for virtual mass was corrected to account for 
the degree of torsional coupling. These data were correlated on the 
basis of four parameters, pitch ratio, diameter, mean width ratio, and 
number of blades. 
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The following conclusions were drawn: 

1. % rsw = 5*775 ^(1.350-PA>) 2 + .228] ^d ^ 3 N* 8 MWR 1 * 13 lbs. (\) 

0.75 P/b 1.25? P/D = Pitch ratio 

d = Propeller diameter (in.) 

N = Number of blades 
MWR = Mean width ratio 

2. Virtual mass as a percentage of total weight of bronze 
propeller blades (blade thickness fraction = 0.05) is as follows: 

P/D 1* Bladed Merchant 3 Bladed DD U Bladed DD 

,9 n% 71# 72# 

1.2 Ul# U2£ 

3. Values of virtual weight under operating conditions were 
significantly less than static test values. 

U. No consistent variation of virtual weight was apparent with 
a variation of J p 

5. Damping factors were about l/3 theoretical values. 



The following recommendations were made: 

1. Compare aluminum and bronze model propellers to investigate 
torsional coupling. 

2. Rim more propellers, including some of varying diameters. 

3. Investigate damping more thoroughly, 

!u Investigate variation of support longitudinal stiffness. 



Thesis Supervisor: Frank Mj Lewis 

Title: Professor of Marine Engineering 
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DEFINITIONS 



Damping constant of propeller and support complete, 
“ a ^a (lb sec) 

T* 7HT 



Damping constant of propeller, B -B ( lb sec ) 

S IftT 

^ , 2 a(W + wt of hub and entrained water) 

Damping constant of support s 

g 

(lb sec) 



Simple theoretical damping constant, ^t (lb sec) 

dJ IftT” 



Torsional coupling factor 
Diameter of propeller, (ft) 
Diameter of propeller, (in) 



Height of bromo-bsnzene and water column (fine scale) for v 
computation, (imi) 

Corrected height of column due to temperature effect, (This is 
negligible in its application to the experiments but the correc- 
tion is a characteristic of the MIT tunnel.) (mm) 

✓v 

Polar inertia, (lb sec~ in) 



v 



Parameter of propeller operation, (dimensionless) 

m , . . . (thrust) / N 

Thrust coefficient, • ( unaensiorf ess ) 



Longitudinal spring constant, 



(lb) 

XinJ 



Torsional spring constant 




(in 11 



Shaft revolutions per minute, (rpm) 
Shaft revolutions per second, (rps) 



Pitch of propeller 3 (ft) 

Pitch of propeller, (in) 

2 

T ^ w 

Slope of plot of W (or weight added) vs, n , a (lb) 

( se ° 2 > 

Observed period of vibration, (sec) 



Natural period of vibration. = 



I 2 2 ' 

\l 1 * g V 
v h r 2 



Temperature of water in tunnel, (°F) 



Water velocity, propeller speed of advance. 0.416^6^17^^ 

(ft) 

(sec) 

Apparent weight of propeller and support complete in water. 

(From period measurement), (lb) 

Apparent weight corrected to condition of infinite torsional 
stiffness of coupling C T W (lb). (W as used in this equation 
was the (average, operating) as tabulated in Tables III - XII) 

Virtual weight of water entrained by blades, W •- W -W..-W (lb) 

6 3. JL S 



Virtual weight corrected to condition of infinite torsional 

l-W 
1 s 



stiffness of coupling, W 1 ~ W -W.,-W_ 



6 

corrected 

Weight of water entrained within hub around shaft, (lb) 



Weight of propeller, (lb) 



Apparent weight of support in water including water entrained 
in housing, nose, and propeller drive collar. (From period 
measurement), (lb) 



Weight added to shaft in excess of nose and collar, 



Reciprocal of time constant, exponential decay. 

(lb sec^) 

(ft 1 *) 



p4 

(sec) 



W h (lb) 



Density of medium. 



INTRODUCTION 



The determination of the effective mass and the damping of marine 
propellers is essential to the analysis of the longitudinal vibrational 
characteristics of a ship’s shafting system. By definition the effective 
mass consists of the actual mass of the propeller plus a virtual mass 
due to entrained water vibrating with the propeller. In the process of 
analysing the shafting system it is necessary to make many aiglneering 
assumptions such as the degree of fixity at the thrust bearing, the amount 
of structure moving with the shaft at the forward end, the yielding of 
the ship’s hull, the inertia of the water entrained by the hull, and the 
effective mass of the propeller. It was our object to contribute to the 
solution of the general problem by making a more accurate determination 
of this effective mass than has been made heretofore. 

The only previous determination was with free-to-rotate nor-driven 
propellers immersed in water and vibrated axially. (I) - Our apparatus 
was designed to determine the virtual mass of propellers under operating 
conditions, that is, rotating and delivering power. The method was such 
that it was also possible to measure propeller damping. 

The basic instrumentation problem was to support the propeller so 
that it was free to vibrate longitudinally (along the propeller axis) 
and so that torque could be applied to drive the propeller. The support 
system as illustrated in Figures II-VT consists of a 12 inch section of 
3/b inch stainless steel shafting supported and longitudinally constrained 
by laminated bakelite bearings at each end of the bearing housing. This 

■^umbers in parenthesis refer to references listed in the bibliography. 
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Figure I 

Operating position adjacent to 
tunnel showing Sanborn Recorder, 
initial-displacement-wire, (arrow) , 
and tunnel window for observation. 





Figure II 

Propeller (Merchant Type) in 
operating position. (Flow is 
left to right. ) Note initial- 
displacement-wire and leads to 
strain gages on the nearest 
rod. 



Figure III 

Exploded view of apparatus 
showing three different 
propellers with associated 
collar (thrust) and nose 
assemblies . 
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Figure IV 

Exploded view of apparatus 
showing collar and nose assemblies, 
thrust collar, torque coupling, 
shaft, support rods, and strain 
gages. 








Figure V 

Apparatus assembled except 
for propeller. 



B'igure VI 

Apparatus assembled and 
installed in propeller 
tunnel (3 bladed DDI4.O9 Type 
Propeller) showing torque 
coupling and initial-displacement 
wire . 
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housing is a section of bronze pipe $ inches long and 1.7 inch outside 
diameter which is itself supported by four 3/16 inch steel rods 12 inches 
long. The rods are bolted at their upper end to 6” x 2" channel beams 
rigidly secured in the propeller tunnel. A collar is keyed to the shaft 
and drives the propeller through dowels. A nose piece is screwed onto 
the shaft to seat against the propeller hub in order to hold the propel- 
ler securely against the collar and to fair the approaching water flow. 

To the driven end of this shaft a torque arm is keyed. This torque arm 
is driven by four springs connected to a similar torque arm on the driv- 
ing shaft. The radius of the torque arm is 2 inches and the effective 
torque radius of the springs is 1*2 inches. Typical tested propellers 
are illustrated in Figure V. 

Propellers having the greatest variation possible in pitch/ diameter 
ratio and mean width ratio were selected for testing. Ore series con- 
sisted of five merchant type propellers, two bronze, three aluminum, hav- 
ing four blades, a MWR of «2ii8, and having P/D ratios from <>8 to 1.2. 

Five aluminum destroyer propellers of three and four blades with MWR of 
,l£S and the same P/D range were selected for the second series. The 
propellers are further described in tables III to XIII in Appendix A. 

To make observations the entire support system was excited by a 
step displacement from the steady state. The signal from, an SR-ii strain 
gage transducer attached to one of the rods was the input to a Sanborn 
Recorder. The recorder produced a trace (see Figures VII and VIII) 
which is a measure of the instantaneous longitudinal position of the 
propeller. The time for each complete oscillation of the trace is the 
damped period of vibration of the system. This period may be simply 
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related to the longitudinal spring constant of the system and to the ef- 
fective mass of the system* By analyzing the variation of the effective 
mass between runs in air and in water it was possible to determine the 
added or virtual mass due to operation in water. The effect of torsional 
coupling was introduced into the results. 

The damping constant is determined by analysis of the exponential 
decay of the displacement trace. Equation 1 was empirically derived to 
describe the curve resulting from a plot of virtual mass versus pitch 
ratio of the propellers. 
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PROCEDURE 



General 

The general approach to the measurement was to compare the natural 
frequency of the vibrating system in water with the natural frequency of 
the system when in air. The frequency and damping measurements were 
determined from the transient response of the system following a step 
displacement from the steady state 0 These measurements were obtained 
from the trace of a "Sanborn Recorder^ Model 150 ,# . The input to the 
recorder was a half -bridge SR-ij. strain gage transducer which was attached 
to one of the vertical rods which constitute the primary longitudinal 
spring system. 

Method of Excitation 

The system was excited by manually displacing the support system 
longitudinally from its steady state position. This was accomplished 
by pulling on a piano wire which was connected to the support system 
through a system of fairleads (Figures I and II). When the piano wire 
was released suddenly^, the system response was an exponentially decaying 
sinusoid o Caution was taken to ensure that after release no tension 
existed on the tripping wire. The longitudinal displacement prior to 
release was approximately l/8 inch. 

Method of Measurement 

Figures VII and VIII show a typical response. The period is deter- 
mined by measurement of the first several cycles of the decaying sinusoid. 
The first cycle was omitted from this measurement. The exact number of 
cycles measured depended upon the rate of damping and upon the "noise 
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level" of the system. This number varied from two to five for runs with 
the propellers installed. The average number of cycles used was four. 

The inverse time constant of the exponential decay (a) was determined by 
use of transparent overlays. Since the damping ratio (a/ w n = aT J 2 if ) 

was always less than ,1 5 no correction was necessary to the damped period 

in order to obtain the natural period. 

The data tabulated for each "Run" represents the arithmetic average 

of the data obtained from no less than five excitations under the same 

operating conditions. The excitations for a given "Run" were made within 

approximately fifteen seconds. 

Calibration of the Support System 

The propeller supporting system was calibrated in air by adding 

known weights to the system in lieu of the propeller. The square of the 

period thus observed (Table I) was then plotted against the weight added 

(Figure IX ). The straight line thus described was then a measure of the 

2 2 

longitudinal stiffness of the spring system. i.e 0 AW/AT Q = S = ljjg/ii if • 

Furthermore the extrapolation of this line back to zero period (Figure IX) 
then measures the weight of the supporting system in air. This weight 
was found to be 1 +. 89 #. 

Each time the system was installed in the tunnel the calibration 
was checked by observing the period of the system with an added weight 
of 6.562#. The ratio of apparent weight to period squared (W /T = S) 

3 l O 

was then calculated (W being ii. 89 # +6.562#). This value of S was then 

€1 

used for the runs made during that day. The slope was checked in the 
same manner at the conclusion of the day's runs. 
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level 11 of the system. Uiis number varied from two to five for runs with 
the propellers installed. The average number of cycles used was four. 

The inverse time constant of the exponential decay (a) was determined by 
use of transparent overlays. Since the damping ratio (a/^ = aT Q / 2 7/* ) 

was always less than .1 5 no correction was necessary to the damped period 
in order to obtain the natural period. 

The data tabulated for each "Run” represents the arithmetic average 
of the data obtained from no less than five excitations under the same 
operating conditions. The excitations for a given "Run” were made within 
approximately fifteen seconds. 

Calibration of the Support System 

The propeller supporting system was calibrated in air by adding 
known weights to the system in lieu of the propeller. The square of the 
period thus observed (Table I) was then plotted against the weight added 
(Figure IX )„ The straight line thus described was then a measure of the 
longitudinal stiffness of the spring system. i.e,AW/AT Q = S =1 . 

Furthermore the extrapolation of this line back to zero period (Figure IX) 
then measures the weight of the supporting system in air. This weight 
was found to be l+.S?#. 

Each time the system was installed in the tunnel the calibration 
was checked by observing the period of the system with an added weight 
of 6. £62#. The ratio of apparent weight to period squared (W /T = S) 

cl O 

was then calculated (W being U<> 8 9 # + 6.£62#). This value of S was then 

2. 

used for the runs made during that day. The slope was checked in the 
same manner at the conclusion of the day's runs. 
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The virtual weight of water and the support damping associated 
with the oscillation of the support system when immersed in water was 
determined by use of a dummy hub equal in size to the hub of the merchant 
type propeller. It was found that for the range of water velocity used 
in the experiments the virtual weight of the support system was indepen- 
dent of water velocity (see Table II and Figure X). The virtual weight 
of the water including that actually contained within the bearing sup- 
port was found to be 0.313# • The virtual weight of the support in water 
with the merchant type nose and collar piece (W ) then equals £.203#. 

This value was used throughout the merchant type tests. The support 
damping was found to be dependent upon water velocity and its relation- 
ship is plotted in Figure XI. 

The virtual weight of water and damping for the support was assumed 
to be the same for the destroyer type propeller installation,, Due to 
changing the nose piece.and collar for the DD type propeller, the virtual 
weight of the support in water (W ) for the DD type propeller is then 

5 

lu*03#. 

Calibration of Recording System 

Since the "Sanborn Recorder" used has paper speed errors as great 
as 1$, it is necessary to correct the measured period for paper speed 
error. The recorder has a time marker which records one second intervals 
on one margin of the recording tape (see Figure VII). Since the timing 
device is an unloaded synchronous motor, the time marker was assumed to 
have only those errors caused by line frequency variations. The line 
frequency was monitered and runs were made only when the line frequency 
was 60cps ± 0.02cps. On this basis corrections were made only for paper 
speed error. 
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Figure XI 

Damping Constant for Support vs. Approach Velocity 
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Method of Calculation 



The value of T which is tabulated in the tables of data was 
o 



corrected for paper speed error only. 

2 



Then : 
and 



W = S x T 
a o 



W e = W a " < M 1 + V 



To correct for coupling between torsional oscillation and longi- 
tudinal oscillation a torsional coupling factor was determined (see Ap- 
pendix B ). 



Then: 



^a corrected ^a 



and 



V = W a corrected * ' W 1 + V 



The value of a which is tabulated on the tables of data was deter- 
mined by means of transparent overlays# 

Then: B ~ 2a VI /g 

a a ° 

and B = B - B 

e a s 



The theoretical damping constant (b) was evaluated from the equa- 
tion (see Appendix B )• 

b = J>nD 3 dK t 
dJ 






RESULTS 



The data are tabulated in Tables III - XIII of the appendix, 
each table giving the data pertinent to and the characteristics of a 
given propeller# Each propeller is identified throughout this paper 
by its respective MIT number# In Figures XII - XV the number next to 
a point indicates the experimental run from which the point was deter- 
mined# 

The torsional stiffness of the system (kg) for the runs above 
200 was 7#3l| ir#/radian and had a ratio of longitudinal frequency to 
torsional frequency of 5> 0 1 for the bronze propellers and It. .It for the 
aluminum propellers# The runs numbered below 200 had a torsional stiff- 
ness k “ 86#5> in#/radian and had a ratio of longitudinal frequency to 

0 

torsional frequency of 1#55> for the bronze propeller and 1#31 for the 
aluminum propellers. The determination of these ratios was only for the 
merchant propellers and is more fully covered in the discussion of re- 
sults section o 

Figure XII gives the determination of virtual mass for the mer- 
chant propellers in the non-rotating static condition plotted against 
P/D ratio# The torsional restraint to the propeller was provided by 
the torsional stiffness of the support and driving system(k^ and by the 
coulomb friction between the system bearings and the system shaft. 

The upper curves in Figures XIV and XV give the resultant virtual 
mass in the non-rotating static condition for destroyer propellers of 
3 and I4 blades respectively# All of these runs were taken with the 
stiff springs of = 86.5> in#/radian 0 
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Figure XIII 

vso p/p^ (Not corrected) 
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Figure XVII 




bsSimple Theoretical Damping Constant 



32 - 





Figure XIII gives the virtual mass in the operating condition 
for the merchant type propellers plotted against P/D ratio. 

The lower curves in Figures XIV and XV give the virtual mass in 
the operating condition for destroyer propellers of three and four 
blades respectively. 

Experimental damping terms (B ) were calculated from the expert- 

© 

mental data and theoretical damping constants (b) were determined from 
the formulation b ^pnd'’ dK^/dJ which is derived in appendix B • Both 
b and B e are tabulated in Tables III - XIII for all propellers where 
the necessary data was available. Curves of the experimentally deter- 
mined B g which correspond to the theoretically determined values b are 
drawn in Figures XVI - XVII for each propeller. 
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DISCUSSION OF RESULTS 



Damping 

The comparison of the theoretical and experimental values of the 
damping constant of the propellers in Figures XVI and XVII shows that 
the experimental values (B ) are about 3C$ of the theoretically 
predicted values (b), The only exception to this was found in propel- 
ler number 81 (bronze) which gave a value of B which was 2 %% greater 
than the predicted value* There was no consistent variation of the 
damping constant with the speed coefficient (J = v^/nd). More theoreti- 
cal analysis is recommended. 

Virtual weight of water without torsional coupling correction 

The virtual weight as determined without the torsional coupling 
correction is tabulated in Tables III-XIII and is plotted in Figures 

XII- XV versus P/D ratio. No consistent variation in W with speed 

© 

coefficient (J) was observed for the operating runs. 

It is to be noted that the value of W determined from the 

e 

operating runs is much lower than the value determined for the runs 
with the propeller stopped. For the merchant type propellers the 
latter values conform very closely to the W as predicted by the Kane 
and McGoldrick formulation (1). 

In general the trend of the values of W as plotted in Figures 

© 

XIII- XV is to decrease with increasing P/D. This is the trend which 
would be anticipated from physical reasoning. If this trend had been 
entirely consistent, it would have been possible to conclude that our 
calculated values of were a true representation of the virtual 
weight. The inconsistency is that the bronze propeller #80 with a 



P/D of 1.195 has a value of W which is appreciably greater than the 
value to be expected from a fair curve through the values of W g deter- 
mined for the aluminum merchant type propellers having lower P/D ratios. 

It is felt that this apparent discrepancy is due to the fact that 
the polar inertia of the aluminum propeller is less than that of the 
bronze propeller and is, because of this lewer inertia, more subject 
to torsional coupling than is the bronze propeller. 

Before investigating the effect of the lower polar inertia of 
the aluminum propellers the relationship between the first approxima- 
tion longitudinal and torsional characteristics of the apparatus were 
investigated. A variation of the torsional characteristics of the sys- 
tem was achieved by varying the stiffness of the torque transmission 
springs. The runs numbered above 200 were run with springs of 

k =7.34 in#/ radian and the other runs were with springs giving a 
© 

k = 86.5 in#/radian. Since this torsional stiffness is much lower 
© 

than that of the drive shaft (k^ = 2500), a first approximation of 
the torsional frequency of the system may be arrived at by considering 
the kg of the torque transmission system alone. The longitudinal 
spring constant has inherently been used to determine the virtual 
mass and is readily evaluated. 

Values of the first approximation torsional frequency and the 
longitudinal frequency are determined below for one arbitrarily selected 
bronze propeller and one of aluminum. As illustrated in Figures XII-XV, 
the shift in the frequency ratios apparently had negligible effect on 
the amount of torsional coupling. (By comparison of 200 runs and other 
runs with the same propeller.) 
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System #2 


System #1 


Huns Above 200 


Runs Below 200 


Soft Springs 


Stiff Springs 


k L = 83.1* #/in, = 7.31* r t d ln 


k x = 89 , k 0 = 86.5 


Bronze Propeller #80 Run 205 


Run 100 


f L = ".lfCS? = 7 ’ h2 /s9 ° 


"L . 1^920 ,#7 ^ //sec 


I ~ (1.25X.0706) = 0.0881; 
P 


I = .0881; 

P 


f Q = 1.15 /sec 


£q - 5.0 /sec 


f L/ f =5.1 
0 


f T / = 1.6 

W f 

0 


Aluminum Propeller #'{h Run 211 


Run 85 


f L = “3533 = 9 ' 66 /sec 


f L = TTSTo = 9 - 9 / 3ec 


(. 28 ? + .25 ). 0706 = 0.038 

t 


~ 7.56 /sec 


i' = 2,21 /sec 

f L/ r - h oh 
0 

Torsional Coupling 


H 

• 

iH 

II 

a> 



To make aHowance for the variation in polar inertia between the 
aluminum and the bronze propellers and to obtain an ’’absolute" value 
of virtual mass, torsional coupling was investigated and a torsional 
coupling factor was derived,, The "absolute" value of virtual mass is 
defined as the value which would have been observed if there had been 



no torsional coupling. 

Torsional coupling is a measure of the interaction of torsional 
and longitudinal vibrations. The greater the torsional coupling the 
less will be the amount of water entrained longitudinally by the rotat- 
ing propeller. 

Applying this to the experiments conducted it may be seen that if 
either the polar moment of inertia of the propeller or the torsional 
stiffness of the driving shaft (or coupling) is low, the amount of 
torsional coupling will be large. 

In the shipboard case (U) little torsional coupling should be 
expected o In our determination the polar moment of inertia and the 
torsional spring constant are low so that some torsional coupling is 
anticipated. Therefore a torsional coupling factor, C^, was derived 
by assuming that the experimental determination corresponded to a case 
of high torsional coupling and that a correction to the condition of 
low torsional coupling was necessary, (i.e. ). 

The basis for this torsional coupling correction is an unpub- 
lished paper by Professor Frank M, Lewis of the Massachusetts Insti- 
tute of Technology. (See extract in Appendix C). 

The derivation of and sample calculations are detailed in 
Appendix B. The values of varied from 1*02 for a bronze propeller 
to 1.11 for an al.umin.um propeller. This correction was applied to the 
system weight W and led to a corrected virtual weight of entrained 

cl 

water W 1 9 recorded in Tables III-XIII and in Figures XVIII-XX. 

Errors in the Measurement 

The major source of error in the measurement was in the assumption 
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that the longitudinal stiffness of the system was constant during the 
runs made on a given day. The stiffness was found to vary by a maximum 
of about 2% from the initial determination made prior to the runs to the 
final determination made after the runs. The reason for this variation 
is unknown. 

The measurement of the period (T ) by the method outlined in the 
Procedure is estimated to be accurate to within 1 %. 

The paper speed correction and the variation line of frequency 
are each estimated to be accurate to within 0,1^. 

The above errors then give an estimated error of 3o2% in the 
determination of the observed period (T )• This corresponds to an 
error of 6,1$ in T and in W • The mean -value of W (aluminum) is 

O cl cl 

about 9# and the mean value of W (bronse) is about 15># for the mer- 

cl 

chant type propellers. The mean value of W for the DD type propellers 
is about 9#, 

This then corresponds to an estimated total error of ,57# iu the 
determination of W (and/or W * ) for the aluminum merchant type propellers 
and for the DD type propellers. The total error in the bronze merchant 
type propellers is then ,96#. This estimate of the error of the measure- 
ment is confirmed very well by Figures XIII-XV, 

Concluding Discussion 

On the basis of the virtual mass corrected for torsional coupling 
an empirical relationship was derived to describe the variation of 
virtual mass with the important parameters. The derivation is given 
in Appendix B and leads to the followings 

W e'(rsw) = 5. 773](l.35 0 - (P/D)f + 0.228~| (^) 3 N 0 * 8 MWR 1,13 (i) 
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Figure XVIII 
average, corrected) vs. P/D 
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.8 .9 1 © 0 1.1 1.2 

P/D 

DDi|09 Type k- Blades 4=11.82" MWR=0. Ij-58 






This equation held for both the low mean width ratio (.2^8) mer- 
chant propellers and the high (,1£8) MWR destroyer propellers. It gives 
the absolute virtual mass, that is the virtual mass with no torsional 

coupling. This corresponds to the shipboard case where torsional 

* 

coupling is not anticipated. The values of W g , VM, ^ e '( rgw )» anc * 
^e(K^M) are ^ abula;ted Table IA. It is our conclusion that 
W 1 / \ gives a quick and sufficiently accurate determination for use 

© ( I* SW j 

in shipboard calculations. 

To improve on our determination, more analysis of the effect of 
torsional spring and inertia coupling is required. The manufacture of 
two propellers, one bronze and one aluminum, of the same characteristics 
would allow comprehensive testing and absolute determination of the vari- 
ation in longitudinal virtual mass due solely to variation in the polar 
inertia of the propeller. Also, propellers of varying diameter should 
be tested to determine the diameter influence and to compare it with 
the theoretical prediction of variation of virtual mass with the cube 
power of the diameter. 



TABLE IA 



Comparison of virtual weight of entrained water as determined 
by various methods. 



W. 



B 



Prop, 

No. 


wt of prop 
corrected to 
Bronze 


P/D 


w 

e 


V V(r-r) 


W e(K£M) 


81 


7.10; 


.796 


2.65 


2.9U 


3.26 


3.3 


72 


7.85 


.90 


1.57 


2.68 


2.63 


3.2 


Ik 


7.85 


.9975 


3*07 


1.99 


2.15 


3.07 


73 


7.65 


1.099 


.87 


1.70 


1.77 


2.95 


80 


7 . h 7 


1.195 


1.83 


2.25 


1.1*7 


2.81* 


30 


10,0 


.890 


2.1*8 


li.31 


1*.18 


7.65 


28 


10.0 


1.11 


1.1|6 


2.72 


2,67 


7.05 


32 


10.0 


1.19 


1.51 


2.50 


2.1a 


6.80 


31 


11.3 


.90 


2.98 


5.30 


5.16 


10.7 


33 


11.3 


1.20 


1.92 


3.26 


3.00 


9.5 




( ^e'(rsw)) 


= 5.775 


[(1.350-p/d) 2 


+ ,228j 


(I5 )3 (« 




W 


= 9100 x — 


1 X 


N x(r~ 


) 3 (MWR) 2 



.23(|) + 1 

where: is the experimentally observed value 

W 1 is the observed value corrected to condition of 
e 

infinite torsional stiffness of coupling 



w i . . 

e (rsw) is the value predicted by equation 

^e (K&l) is the value predicted by the Kane and 
McGoldrick formulation (l) 



CONCLUSIONS 



1. Virtual weight of entrained water can be expressed by the 
following equation: 



V 

(rsw) 



= 5.775 |jl.350-P/D) 2 + .22^ lbs. 

0.75 < P/D < 1.25 



P/D = Pitch ratio 

d = Propeller diameter (in) 

N = Number of blades 
MWR = Mean width ratio 

2. The following percentages represent the virtual weight for a 
sampling of bronze propellers on the basis of total weight of propeller 
blades: (blade thickness fraction — 0.0f>) 

P/D U B laded Merchant 3 B laded DP It Bladed DP 
.9 71 # 73 # 12 % 

1.2 itlt# It# h2% 



3, The values of virtual weight under actual operating conditions 
were significantly less than the values predicted on the basis of non- 
rotating static tests. This was particularly true for the DD types* 

It. There was no consistent variation of virtual weight with varia- 
tion of J in the operating range, 

5>. Magnitudes of damping factors under operating conditions were 
about 30 % of the theoretically predicted values. There were no signi- 
ficant trends with a variation of J or v . 
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RECOMMENDATIONS 



1, Investigate torsional coupling factor to a greater degree, 
particularly the effect of varying polar moment of inertia. It would 
be advantageous to compare aluminum and bronze propellers having the 
same geometrical characteristics, 

2, Obtain data for more propellers to establish curves more 
thoroughly. To determine or prove scaling factor run geometrically 
similar propellers of varying diameters, 

3, Investigate damping more thoroughly varying the parameters, 
lw Investigate the variation of longitudinal stiffness of the 



support 



APPENDIX 



A. Summary of Data and Calculations 

TABLE l 



CALIBRATION OF SUPPORT SYSTEM IN AIR 



RUN 


Wt 


T o 


T 




added 


0 


o 


1 


6.^62 


.111*2 


.01320 


2 


2.1*69 


.0919 


.0081*1* 


3 


3.1*69 


.0993 


.00986 


U 


0.375 


.0775 


.00600 



The weight of the supporting system in air including the nose and 

collar used with the merchant type propellers was determined to be 1 *. 89 #, 

This was determined by extrapolating the curve of weight added versus 
2 2 

T q to the value of T q equals zero. (See figure IX ). 
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TABLE II 



RUN N 



1 0 

2 0 

3 0 

it 0 

5 0 

6 0 

7 0 



CALIBRATION OF SUPPORT SYSTEM IN WATER 



Added weights Dummy hub wt = .812# 







Water 


in hub = , 


.375# 






Wt added 


= 1.187# 


h 


V 


2 

T T 


W 


W 


ca B 




0 


o 0 




s 


s 


0 


0 


.0883S .007795 


6.1*3 


5.21t3 


1.5 .595 


32.9 


2.39 


, 087 it 0 . 00761*0 


6.30 


5.113 


2.0 .799 


35.9 


2. £0 


.08795 .007735 


6.38 


5.193 


2.0 .79it 


ll.it 


1.1*1 


.08820 .007780 


6 . 1*2 


5.233 


1.7 .675 


17.0 


1.72 


.O 88 OS .007750 


6 . 1*0 


5.231 


2.0 .791* 


1 * 1.6 


2.69 


.08735 .007630 


6.30 


5.221 


2.0 .791* 


66 .it 


3.1*0 


.08725 .007610 


6.28 


5.219 


2.1 .831* 



W = W-1.187# 

S 
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TABLE III 



M.I.T. PROPELLER NO. 81 



Merchant Type 




w = 7.1*37# 


Material 


= Bronze 


O 

1 ! 

5^ 


Diameter (d) 


= 11 . 90 ” 


W, = 7.1*37# 


No. of Blades (N) 


= 1 * 


JL 

w,, = 5.203# 


Pitch Ratio (P/D) 


= 0.796 


S 



Mean Width Ratio (MWR) = 0.21*8 



dK t /dJ .!*< J <.9 = .1*52 



= 836 (Runs 91 “ 9h); 86 (Runs 215 - 217) 



RUN 


h 


t 


V 


N 


j 


T 


T 


W 


w 


a 


B 


B 


B 


b 






0 






0 


0 


a 


e 


a 


s 


e 




9 1 


0 




0 


0 


mm 


.1380 


.01903 


6.910 


3.27 


0.8 


.99 


69 


0.1*0 




9 2 


17.0 


696 


1.72 


61 


.68 


.1361* 


.01861 


6.560 


2.92 


3.0 


2.85 


.71* 


2.11 


2.6 


93 


336 


696 


2.91 


62 


.95 


.1366 


.01865 


6.590 


2.95 


1*.0 


3.80 


.79 


3.01 


2.17 


9 k 


60.0 


696 


363 


160 


1.21 


.1327 


.01761 


6.730 


2.09 


1*6 


1*.28 


.82 


3.1*6 


2.28 


215 


0 


- 


0 


0 


- 


.1361* 


.0186 


6.6 


261 












216 


316 


70 


2.3U 


63 


.98 


.1333 


.01775 


6.50 


1*86 












217 


25 


70 


2.08 


61* 


.87 


.1333 


.01775 


6.50 


1.86 













Run 217 readings of T 0 were erratic. 



W (average, operating) = 6.976# 

cl 



W ' = 2.9lt# 
e 
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TABLE IV 



M.I.T. PROPELLER NO. 72 
Merchant Type 
Material 
Diameter (d) 

No. of Blades (N) 
Pitch Ratio (P/D) 

Mean Width Ratio (MWR) 

dK t /dJ .5<J<1.1 
S 



RUN 


h 


t 


V 


N 


J 


T 


T 


w 


W 


a 


B 


B 


B 


b 








0 






0 


0 


a 


e 


a 


3 


e 




87 


0 




0 


0 


— 


.1137 


.01292 


10.80 


3.281; 


1^6 


1.07 


.59 


0 .U 8 




88 


18.6 


69.5 1.79 


151 


.71 


.1022 


.OlOU; 


8.73 


1.211; 


1.5 


0.85 


.75 


0.10 


2.21 


89 


30.3 


69.5 


; 2.29 


15>2 


.91 


.1050 


.01103 


9.22 


1.70U 


3.5 


1.98 


.78 


1.20 


2.23 


90 


59.0 


69.5 


> 3.20 


1^6 


1.23 


.1056 


.01115 


9.32 


1.90U 


2.5 


1.1a 


.82 


0.59 


2.29 


213 


0 


- 


0 


0 


- 


.1111 


.0123U 


10.06 


2.UU 












21 h 


Uo 


70 


2.61; 


ll ;0 


1.13 


.1060 


.on23 


9.16 


1.61; 













W (average, operating) = 9.107# 

3 . 

W ' = 2.68# 

6 



Aluminum 

11.90" 

h 

0.900 

0.2U8 

•U65 



W = 2.250# 
P 

W h = 0.063# 



W ]L = 2.313# 
W s = 5.203# 



= 836 (Runs 87 - 90 ); 815 (Runs 213, 211;) 
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TABLE V 



M.I.T. PROPELLER NO. lh 





Merchant Type 






w p = 2.250# 








Material 




= Aluminum. 


w h = 0.125# 








Diameter (d) 




= 11.90” 


















W. = 2.375# 








No. of Blades (N) 


= h 


-L 
















W, = 5.203# 








Pitch Ratio (P/D) 


= 0.9975 


5 








Mean Width Ratio (MWR) 


= O.2U8 










dK t /dJ 


.5 <1.2 


CM 

r- 

-=t 

• 

11 










S 






= 836 (Runs 83 - 


86); 815 (Runs 210 - 212) 






RUN 


h t 


V 


N 


2 

J T T 


W W „ B B 


B 


b 






0 




0 0 


a e a a s 


e 




83 


0 


0 


0 


- .1133 .01281; 


10.73 3.152 1.8 1.20 0.59 


0.61 




8U 


32.0 68. 


0 2.36 


150 


, 9 h .lOll .OIO83 


9.05 1.U72 3.1 1.66 0.78 


0.88 


2.23 


85 


57.3 68. 


0 3.16 


152.5 1.2U .1010 .0102 


8.53 .952 2.U 1.29 0.82 


0.1*7 2.27 


86 


16.9 68.0 1.711; 


150 


.68 .1001 .01002 


8.38 .802 1.2 0.6U 0.73 


0.09 2.23 


210 


0 


0 


0 


.1063 .01130 


9.21 I.63 






211 


36 ?0 


2.50 


1U8 


1.01 .1035 .010712 


8.71 1.16 






212 


1*7.5 70 


2.88 


UU2 


1.21 .1031 .010630 


8.67 1.09 







W (average, operating) = 8.67!*# 

Si 

vf e ' = 1.99# 
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TABLE VI 



M.I.T. PROPELLER NO. 73 



Merchant Type 
Material 
Diameter (d) 

No. of Blades (N) 

Pitch Ratio (P/D) 

Mean Width Ratio (MWR) 



= Aluminum 

= 11.90" 

= 1 + 

= 1.099 

= 0 . 21+8 



W p = 2.219# 
W h = 0.125# 

W 1 = 2.33U# 
W s = 5.203# 



dK t /dJ ,6£<J*1.2 = .1+72 





z 

S 








— 


881 (Runs 102 


- 105)3 815 (Runs 208, 


209) 


RUN 


h 


t 


V 


N 


J 


T T 2 


w w 


a 


B 


B 


B b 








0 






0 0 


a e 


a 


s 


e 


102 


0 




0 


0 


__ 


.1071* .01151* 


10.16 2.623 


1.7 


1.07 


.59 


.1*8 


103 


19.5 67.5 


1.81+ 


151* 


.72 


.09706 ,0091*2 


8.30 0.773 


2.3 


1.17 


.75 


.1*2 2.29 


10l+ 


38.0 67.5 


2.57 


156 


.99 


.09668 .00935 


8.21* 0.713 


2.7 


1.1*0 


.80 


.60 2.32 


105 


61.0 6?.5 


3.26 


152 1.28 


.09938 ,00988 


8.70 1.173 


2.2 


1.15 


.83 


.32 2.26 


208 


0 


mm 


0 


0 


- 


.1085 .01177 


9.59 2.05 










209 


1*5 


68 


2.79 


150 


1.12 


.101*0 .01082 


8.82 1.28 











W (average, operating) = 8.pl5# 

9 . 

W e ' = l.?o# 
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TABLE VII 



M.I.T. PROPELLER NO. 80 



dK t /dJ •7<'J<1.3 



Merchant Type 




W = 7 . 1 * 7 C# 

.H 


Material 


= Bronze 


0 

it 


Diameter (d) 


= 11.90” 


w, = 7 . 1 * 70 # 


No. of Blades (N) 


= u 


± 

= 5.203# 


Pitch Ratio (P/D) 


= 1.195 


0 



Mean Width Ratio (MWR) = 0.21*8 

= .1*91* 

= 881 (Runs 96 - 100)5 8l5 (Runs 20l* - 207) 



RUN 


h 


t 


V 

0 


N 


J 


96 


0 




0 


0 




97 


0 


- 


0 


0 


- 


98 


16.0 


67 


1.77 


150 


.708 


99 


35.0 


67 


2.1*6 


150 


.985 


100 


57.5 


6? 


3 .16 


150 


1.261* 


201* 


0 


- 


0 


0 


- 


205 


59 


68 


3.20 


151* 


1.25 


206 


57 


68 


3.15 


152 


1.21* 


207 


1*1 


68 


2.67 


105 


1.53 



W. 



W 



a 



B * B o 

a s 



B 



e 



- 01761 1S - 51 2 Mo .8 0~77 .59 7l8 

.12728 .01621 1U.30 1.627 2.3 2.07 .71* 1.33 2.33 

.12828 . 0161*6 ll*.5l I.837 2.1* 2.12 .79 1.33 2.33 

.12920 .01669 U*.71 2.01*7 3.1* 3.10 .82 2.28 2.33 

.1375 .01890 1$. 1*1 2. 71* 

.13l*8S .01818 ll*.8l 2.11* 

.1352 .01828 lii.90 2.23 

.131*3 .01805 1U.71 2.01* 



W & (average, operating) = — g— = ll*.656# 
W e * * 2.25# 
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TABLE VIII 



M.I.T. PROPELLER NO. 30 



DDl*09 



Material 
Diameter (d) 

No. of Blades (N) 
Pitch Ratio (P/D) 



= Aluminum 
= 11 . 82 " 

= 3 

= 0.890 



Mean Width Ratio (MWR) = 0.1*58 



dK t /dJ .!*< J < .93 = .U27 

S = 881 



W = 2.875# 

JT 

W h = 0.063# 

W 1 = 2.938# 
w. = 1*.5 o3# 

0 



RUN h t 



v 

o 



N J 





W 

e 





B 

e 



b 



106 

107 

108 

109 



0 - 0 0 .12136 .01172 13.00 5.56 

18.0 67.5 1.77 152 .70 .10630 .01130 9.96 2.52 

37.5 67.5 2.55 150 1.02 .10610 .01126 9»9h 2.50 

60.0 67.5 3.23 150 1.29 .10590 .01122 9.90 2.1*6 



1.52 1.23 .59 
1.55 .95 . 71* 
2.1*2 1.1*9 .79 
2.1*5 1.51 .82 



•61* 

.21 2.00 
.70 1.98 
.69 1.98 



(average, operating) = 9.933# 

cl 



W e ' = li.31# 
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TABLE IX 



3 Aluminum 

3 11 * 82 ” 

= 3 

3 1.11 
= 0.1*58 
= .1*81* 

3 881 



W p 3 2.875# 
W h 3 0.12$# 

W 1 3 3.000# 
W g = ii.503# 



M.I.T. PROPELLER NO. 28 

DD 1*09 

Material 

Diameter (d) 

No. of Blades (N) 

Pitch Ratio (P/D) 

Mean Width Ratio (MWR) 

dK t /dJ .65 <J <1.15 

S 

RUN h t v N J T 
o o 



111* 


0 


- 


0 


0 


- 


.1166 


115 


17.5 


68 


1.7 U 


150 


.70 


.10075 


116 


39.0 


68 


2.60 


152 


1.03 


* 


117 


57.0 


68 


3.1 U 


15 H 


1.22 


.1022 


118 


1*0.0 


68 


1.63 


150 


.65 


■* 



T 


w 


W 




B 


B 


0 


a 


e 


CL 


a 


s 


01359 


11.97 


1*.1*7 


3.75 


2.78 


.59 


01015 


8.96 


1.1*6 


1.70 


.96 


.71* 








6.62 


3.72 


.80 


oiol * l * 


9.20 


1.70 


2.5 U 


1.1*5 


.82 








5.50 


3.09 


.73 



* Damping is too severe for period measurement. 



W„ (average, operating) 3 9.080# 



w / = 2.72# 



B b 
e 



.19 

.22 2 . 21 * 
.92 2.27 
.63 2.30 
.36 2.21 



TABLE X 



Mol .To PROPELLER NO. 32 



DD 1*09 

Material = Aluminum 



Diameter (d) =11,82" 

No, of Blades (N) =3 

Pitch Ratio (P/D) = 1.19 

Mean Width Ratio (MWR) = 0.1*58 
dK t /dJ .7<J<1.2£ = .1*1*1 

S = 881 



W p = 2.78c# 

w h = 0.063# 

= 2.81*3# 
W g = 1*.503# 



RUN h t v N J T 
o o 




a B B B b 
a s e 



123 

12l* 

125 

126 



0 - 0 0 .10650 
18.5 68 1.79 U*9 .72 .0991*8 
38.0 68 2.57 150 1.03 .10000 
57.1 68 3.15 151 1.26 .10105 



0113li 10.00 2.651* 3.1*7 2.16 

00989 8.70 1.35 1.80 .99 

01000 8.81 1.1*6 2.50 1.37 

01021 9.00 1.65 3.06 1.68 



•£9 1.57 
.71* .25 2.03 
.80 .57 2.01* 
.82 .86 2.06 



(average, operating) = 8.837# 

d 



w ' = 2.50# 

6 
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TABLE XI 



M.I.To PROP E LLER NO* 31 
DD U09 
Material 
Diameter (d) 

No* of Blades (N) 
Pitch Ratio (P/D) 



= Aluminum 

= 11 . 82 " 

= U 
= 0.90 



Mean Width Ratio (MWR) = O.U58 



S 



= 881 



W p = 3.25C^ 
w h = .062# 

w x = 3.312# 
W = U.503# 

O 



RUN 


h 


t 


V 


N 


j 


T 


T 2 




w 


rr 


B 


B 


B 








0 






0 


0 


a 


e 


U» 


a 


3 


e 


no 


0 




0 


0 




•129U 


.0167S 


lU.76 


6.9US 


1.90 


1.711 


.59 


1.15 


111 


i5.o 


67.5 


1.61 


15U 


.63 


.1090 


.01187 


10. US 


2.63 


1.80 


1.21 


.73 


•U8 


112 


36.0 


67 .S 


2. $0 


1S0 : 


1,00 


.no7 


.01226 


10.79 


2.97 


2.28 


1.53 


.79 


.7U 


113 


58.5 


67.5 


3.18 


1S3 : 


1.2s 


.n23 


.01261 


n.n 


3.29 


2.30 


1.5U 


.82 


.72 



K, curves not available. 



W (average, operating) = 10.783# 

a 

w e J = 5.30# 
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TABLE XII 



M.I.T. PROPELLER NO. 33 



DD U09 




Wp = 3.2SC# 


Material 


= Aluminum 


w h = 0 . 063 # 


Diameter (d) 

No. of Blades (N) 


» 11.82" 


W T = 3 . 313 # 


= U 


J- 






w, = u.503# 


Pitch Ratio (P/D) 


= 1.20 


9 


Mean Width Ratio (MWR) 


= 0.1£8 




S 


= 881 





RUN 


htvNJT T 2 


W W 


_ B 


B 


B 




1 

! 

1 o 
1 

1 

1 O 

o 


a e 


a a 


3 


e 


119 


0 - 0 0 .1275 .01626 


Uu31 6.U9U iwl2 3.67 


.59 3.08 


120 


16.0 68 1.67 150 .67 .10565 .011162 


9.8? 2.03 


1.78 1.07 


.73 




121 


3U.5 68 2.U5 152 .97 .10369 .010751 


9.1i9 1.67 


2.0? 1.2l» 


.79 


.15 


122 


62.0 68 3.28 ]£l 1.31 .10590 .011218 


9.90 2.08 


2.il8 1.50 


.83 


.67 



Kj. curves not available. 

W (average, operating) * 9.7V?# 

3 * 

w e ' = 3.26# 



TABLE XIII 



M.I.T. PROPELLER NO. 35 



RUN 



DD JU09 




W p = 0.750# 


Material 


= Aluminum 


W h « 0.030# 


Diameter (d) 


= 8.001" 








w n = 0.780# 


No. of Blades (N) 


= 3 


JL 






W = lt.lii3# 


Pitch Ratio (P/D) 


= 1.07 


s 


Mean Width Ratio (MWR) 


= O.U58 




S 


= 907 




h tv N J 


T T 2 W 


W a 


0 


0 0 a 


e u 



130 1 0 - 0 0 .08U79 ,00720 6.5h 1.617 1.3U 

130 17.0 68 1.72 232 . hh . 0669 ^ 

131 18.5 68 1.79 150 .72 * 

132 20.0 68 1.86 328 , 3 h .07022 

133 1H+.5 68 iui;6 330 .81 .06975 

13k Ui.O 68 2.76 330 .50 .06896 

135 11.0 68 .U; 330 .08 .0712^ 



k .00li90 U.15 '.h7 




■^Damping is too severe for period measurement. 



A negative W is not possible so this set is not included for 
analysis. The thrust collars were not properly seated against the 
bearings. 
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B. Sample Calculations 



Torsional Coupling Factor 
Derivation 

1. Referring to Appendix Professor Lewis' unpublished notes, 
equation , we see that if we neglect the damping terms we arrive at 
the following: 



assume 

where 



("l* m e • m e 2) * “ 2 K' k e + 1c l V + h k e = 0 

\"9 = (“l' m 9 - m c 2) 

m^ = W a /g = measured apparent mass 



, _ ^[a (corrected) _ 
g 



m^' = ~ ' ' = apparent mass corrected for torsion- 



m 



al coupling 
fl = I + pi 1 (effective polar inertia) 

y p p 

j3 ^ 0.2$ to include torsional virtual inertia 



I = polar inertia of propeller being tested 
B 

Ip* = polar inertia of bronze propeller of same dimensions 
as actual propeller 



m 



= W, 



2 . 



p = pitch, in. 

Suppose k Q * 0 (close to our k Q = 7.3U 

“ U K’ m e * m c 2) ‘ “ 2 ra e \ = 0 




m, 



e 
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3. Suppose kg = OO 




lu Correction of frequencies to k 0 = 



Define = torsional coupling factor 

V = V =c_ 

2 — I 

<*2 \'-™ c 




Sample Calculations 

1. Bronze Propeller #80 P/D - 1.2 

Ip = nD^bt # sec^ in (Reference 6) 

where n = no. of blades 

D = diameter, in 

b = maximum developed blade width, in 

t = maximum developed blade thickness at 1/2 radius, axis to 
tip 



_ MWR D 

“7BH2 



(Reference £ 
page 1*>7) 



MWR 

b 




”L 



Mean Width Ratio 

= 3<S1 ^ 



0.001*6 

“385 



g o 



(U)(11.90) 3 (3.£L)(jj) = 0.0706 



lk.656 

“1BT“ 



0.0380 



m Q a (1 + 0)1 = (1.25) (.0706) = 0.0881* 
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m 



- »= 2 -m^ --"nr- ^ 



V = "t (1 + EE“) = -t c i = \ (1-01803) 



Ve 

Cj = 1*01803 

2. Aluminum Propeller #72 P/D = 0*9 

W 



"L = i7 = “53S 



m_ = I + 

0 p 



= = °- 02 36 

f*'p = ( S oi bronze ) V * **, 
* (0.287 + *25) (.0706) = 0.0379 



rn 



= = 0.0103lt 



Torsional Coupling Factor 
m^* = m^ (1 + m c 2 



"L m 0 



) = “t 



Cj. = 1.1197 
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Determination of Empirical Equation 



The term for virtual mass, W 1 has four major parameters as follows: 

0 

d (diameter) (inches) 

N (number of blades) 

MWR (mean width ratio) 

P/D (pitch ratio) 

In equation form W * can be expressed as follows: 

6 



V = C rsw f(d) f(N) 



where C is a constant. 
/ rsw 



Analysis of figure by trial and error for P/D term. 



7 



Z ^ C = y/x^ where 

X 



r _ w • 

\y - e 

\ vH 

1 A 



y c= .85 x- 


X 


c 




= 2.8 


X 


2 

X 


Z 


.1 


.01 


.028 


.2 


♦Oil 


.112 


.3 


ON 

o 

• 


.252 


•h 


.16 


CO 

• 




.25 


.700 



vn* 

e P/D=l 
x = 1.350-P/D 



- 0.600 



• » 



Shift this curve up O.OU. 

- 0.6it0 



r - w * 

7 = e 
VT 



6 P/D=l 
X a 1.350-P/D 
This best matches plotted values. 



y = 2.80 x 
W * 



e 

VM 

e 



- 0.6ii0 « 2.80 (1.350-P/D)' 



P/D=l 



V 

w 



P/D=l 



2.80 (1.350-P/D)^ + 0.6U0 

.75 < P/D <1.25 



( 2 ) 
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The first term f(P/D) is of the form of equation (2)* Dividing 



by 2*80 simplifies the expression, 

W * = C ftl^O-P/D) 2 + 0.228 f(d) f (N) f(MMR) 
e vsw l J 

The term f(d) is presumed to be (d^) since W * should vary as the volume* 
To simplify this expression divide d, the diameter of the propeller in 
inches, by 10* 

V = C rsw ft 1 * 350-P/D ) 2 + 0,22oj ^) 3 J f(N) f()WR) 

The term f(N) is presumed to be in the farm of N a . This is determined 
from the following comparisons* 

W * 

vJ-r^ * = (y) a = 1.225 a = 0.705 

0 30 



W « 

e 



5 rr 1 = rfs = ( 5 )4 = 1 * 3 ° 5 a = °- 92S 

e 32 



Select a = 0*8 



W * = C 
e rsw 



j(1.350-P/D ) 2 + 0.220] J(^) 3 ] ^(N)° 4 j f(nm) 



The term f(MWR) is presumed to be in the form of (MWR) * This is deter- 
mined from the following comparisons: 



W » 



e 2 g (corrected to 1* blades) = (2.72)(ty) = 3.1*2 



1*P 8 _ 



W « 

e 



pr^ 1 = £73 = = 2 *oi° * - 1.11,0 

e ?3 



W 



W* 



31 _ 5.30 _ ,.l£8a b _ . _ - ,. n 

"* b - l*9o5 b » 1.120 



W« 72 2^7 



Select b = 1,13 
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1.0 1.1 1.2 1.3 







(1) 



V = C rsw [(1.350-P/D) 2 + 0.228^ [(yg) 3 ] [o0° ,8 j [(MWR) 1,13 

The value of the constant is obtained by substituting the data into 
equation (1)., 

Propeller 81 

V = 0 rro <*S3l*)(1.686)(3.03)(.2'07) = C psw 

Propeller 72 

W e' = 0 rsw (.li31)(1.686)(3.03)(.2l57) = J.5S C rsIf 
Propeller 7k 

VM C rsw (.353)(1.686)(3.03)(.2Crr) = .373 C raw 
Propeller 73 

w e * » C rgw (.290)(1.686)(3,03)(.207) = .306 C rgw 
Propeller 80 — W * is obviously too high. 

Propeller 30 

W e' “ 0 r ^,(.U3?)(1.6Sl)(2.Ul)(J£ll.) = .m c raw 
Propeller 28 

W e' = C rro (*281)(1.651)(2.1a)(4ial.) = Mi c rglr 
Propeller 32 

M e ' = C r3M (.2536)(1.651)(2.U)U01l) = .1(18 
Propeller 31 

W e ' = °r aw (. 1 t31)(1.6Sl)(3.°3)(Jal() =■ Mb 0 rOT 
Propeller 33 

V 35 c r 3 W (* 2 ^)( 1 * 6 ^ 1 )(3.o3)(Ja)i) = .520 c rsw 
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Prop. 

No. 

81 

72 

7l*. 

73 

00 

30 

28 

32 

31 

33 

( 1 ) 



W ' 

e 


W * 

e 

5 — 

raw 


^raw 


W ' 
e 


2.9U 


.565 


5.19 


3.26 


2.68 


.1*55 


5.87 


2.63 


1.99 


.373 


5.33 


2.15 


1.70 

2.2£ 


.306 


5.56 


1.77 


U.31 


.723 


5.96 


1*«16 


2.72 


.1*63 


5.87 


2.67 


2.50 


.1*18 


5.98 


2.U1 


5.30 


.8?1* 


5.93 


5.16 


3.26 


.520 


6.27 


3.00 




0 

raw 

ro AT 


- ?1«*6 
9 ' 


- 5.' 



raw 



5 . 775 




5.773 £(1.350-P/b) 2 * 0.228] [(j§-) 3 J J(N) 0#f jj j(HWR) 1#1 ^ 
•75o<p/fo *1.25 
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Derivation of Simple Damping Constant 



b = damping constant ( - • 



= thrust coefficient 
thrust (lbf) 



t 
T 

vP = density = 1.938 for F.W. 

ft 4 

n = revolutions per second 
dT 



b = 



dv 



_ T 



J>n 2 O h 


1 


V 


= o 
nD 


V 

0 

r\ 1 


_ dT_ 


-j> n st. 


dv 

o 


* nD 



2jU t 



OJW 

dT~ = i/ ,nD 



dK, 

i 

dJ 



In the above q) (partial derivative of) is indicated by d. 
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C. EXTRACT FROM UNPUBLISHED PAPER 



THE COUPLING OF LONGITUDINAL AND TORSIONAL VIBRATION 
OF MARINE PROPELLER SHAFTING 

Unpublished paper by 23 February 195U 

Frank M, Lewis, Professor, Massachusetts Institute of Technology 

The longitudinal and torsional oscillations of marine propulsion 
systems are ordinarily calculated on the assumption that they are inde- 
pendent phenomena. In actuality they are coupled through the intermedi- 
ary of the water reactions, so that a torsional excitation produces both 
longitudinal and torsional motion, and longitudinal excitation does like- 
wise. 

If the longitudinal and torsional frequencies are well removed from 
each other this coupling effect will be of only small importance, but if 
a longitudinal and a torsional frequency are close or coincide, the effect 
is very important altering both the position and amplitude of the criticals. 
Anticipating the results to a certain extent it is noted that for coinci- 
dence of the frequencies the virtual inertia and the hydrodynamic damping 
vanish for both longitudinal and torsional vibration. That this should 
be so is obvious on consideration. Assume a propeller consisting of heli- 
coidal surfaces of zero thickness. For equal longitudinal and torsional 
frequencies a possible mode of vibration consists of a coupled torsional 
and longitudinal motion in which the propeller screws itself back and forth 
through the water. There would be zero virtual inertia for such a motion 
and the damping would be limited to surface friction, a very small effect. 
While for a propeller of finite thickness the inertia and damping effects 
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would not actually vanish, they would become small. 

Passing to the derivation of the general relationships it is noted 
that there are two types of coupling, virtual inertia coupling and damp- 
ing coupling. 

The Virtual Inertia Coupling 

There are three virtual inertia factors, these are j 
T he virtual inertia for pure torsional motion, the virtual inertia 
for pure longitudinal motion and the coupling virtual inertia. These are 
not independent factors and an approximate relationship between them can 
be derived as follows :• — 

Assume a propeller consisting of helicoidal surfaces. In Fig. (1) 
let A B be a blade element at angle to the plane of the propeller. 




Let this be given a longitudinal displacement and a tangential dis- 
placement r 0 o Resolving these into components parallel and perpendicular 
to the blade elements, there is found to be a normal displacement 
n = x cos OC - r 0 sin OL 
and a parallel 

p = x sin (X + r 0 cos OC 

Only the normal motion will have virtual inertia attached to it so that 

the normal force can be written: 

dF = (dm + dm* )n 
n 

and the parallel force 
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= drop 



where dm and dm' are the inertia mass in air and the virtual water mass 

associated with a blade element® Further resolving F and F into forces 

n p 

parallel and perpendicular to the axis there is found. 

dF = dF cos d+ dF sin 0( 
x n p 

dQ = rdF„ = rdF cos 0(- rdF sinCX 
0 p n 

substituting and reducing there is obtained 

dF = dirK + dm* (cos^CXx - sin C/cosO(r*0*) 

dQ = dmr^0 + dm' (sin^O(r^ 0 - r sinO<cosC<*x) 

Summing overall elements* 

F x = rax + x ^ dm* cos ‘^0( - *0 ^dm' r cosO(sinO( 

* • 2 2 

q =: „ x ^"dm* r s±aO(cosO(+ 0 ^r sin CXdm' 

where m and I are the mass and polar inertia constants W/g, J/g of the 
propeller in air, but excluding the hub* 

♦ • 

Considering the torque expression we note that if x = 0 there is 
pure torsional motion* It is known then that the virtual inertia is 
equivalent to an increase of the polar inertia of the propeller of the 
order of 25$. Let this fraction be J3. This method of expressing the 
virtual inertia as a fraction of the propeller inertia in air is crude 
but will suffice for the moment* 



There is thus obtained 

dm’ r^ sin^ = pi 

Still considering the torque equation it is noted that the virtual 
inertia terms must vanish if x and 0* are so related that the propeller 
moves parallel to its helicoidal surface. This condition is expressed by 
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where P is the pitch. Equating the dm 1 terms to zero with this relation- 
ship there is obtained 



considering the F expression the dm'* terms must likewise vanish if 

H Jv. 



This leads to 

J^dm 1 cos ^0( = ^^^dm* r sin OC oos <X = 

This should be the virtual inertia for a pure longitudinal motion. Inser- 
tion of numerical values shors that it comes very close to the value, 
in itself a very crude approximation, which is ordinarily assumed. The 
expressions for the longitudinal force and the torque now become 



These are the final expressions for the inertia! .forces associated 
with the propeller. 

For pure longitudinal and torsional motion these reduce to the values 
ordinarily used and the virtual inertia component becomes zero for motion 
parallel to the helicoidal surface. 

The expressions can be more conveniently written 




P 





F = m~ x - m 9 
x Ti c 



Q = m. 9 - m x 
9 c 
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when m^, m^, m c are the longitudinal, torsional, and coupling inertia 
coefficients respectively. 

It is now assumed that there exists a simple harmonic longitudinal 

motion x of frequency « superposed upon a rotary simple harmonio motion 

0 of frequency w. Exciting forces and torques F and 0 e of frequency w 

0 



constants and neglecting damping the veotor farce and torque equations 
for a single mass system aret — 



The 2 natural frequencies of the system are given by the roots 'of 




F - k T x+ itl. 



e L 1. 




2 2 

Q - k Q 0 + iru w 9-m ^ x = 0 

e©0 c 

These equations have the solution 



where 




the above equation, i*e. A 3 0. 



With the substitutions 




P 



m 9 = I (1 + p) 




The frequency equation reduces to 




r 



p 
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Considering the A = 0 equation as applied to a single degree of 
freedom system it is evident that in the absence of the coupling term 



nal and torsional frequencies with water inertia included. One of these 
will be higher than the other. The effect of the coupling inertia m Q is 
to lower the lowest of the two frequencies and raise the higher. It 
will be usual in marine installations for the lowest torsional mode to 
occur at a much lower frequency than that of the lowest longitudinal 
mode, 

(Author's Note: 

The above is a liberal extract from the original paper by Professor 
Frank M. Lewis. It is hoped that there has been no distortion of the con- 
tents by using for the most part only those portions of the paper which 
dealt with our problem.) 



the roots would be w. 



i = ~ 

L ^ 



and i.e., the separate longitudi- 
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APPENDIX D 



Instrumentation Data 
Recorder* Sanborn l£0 

Model 1^2-UOO Serial number 18. 

Gages. SRU, Baldwin Lima Hamilton 

Type A- 7 

Resistance 120.0 ± 0.3 ohms 
Gage Factor 1.91 ±2% 

Lot No. 232 - 11B72 
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